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Abstract 

Coupling between electronic and mechanical degrees of freedom in a 
single electron shuttle system can cause a mechanical instability leading 
to shuttle transport of electrons between external leads. We predict that 
the resulting low frequency current noise can be enhanced due to ampli- 
tude fluctuations of the shuttle oscillations. Moreover, at the onset of 
mechanical instability a pronounced peak in the low frequency noise is 
expected. 
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1 Introduction 



Charge transport in nanostructures is a major research area, both theoretically 
and experimentally. Apart from the average current flowing through a structure 
in response to an applied external field, fluctuations and correlations in time 
of this current are of interest. By studying the current noise power spectral 
density (PSD), information about the charge transport process can be extracted 
that may not be accessible through studies of the average current alone. P An 
example of this is shot noise, arising due to the discreteness of the charge carriers, 
(electrons). [2] 

In Nanoelectromechanical Systems (NEMS) (3j|, mechanical degrees of free- 
dom affect and/or are affected by charge transport through the device. One 
such system is the single electron shuttle system 0], also known as the nanoelec- 
tromechanical single electron transistor (NEM-SET). The system consists of a 
metallic grain embedded in an elastic material between two bulk leads, forming a 
mechanically soft Coulomb blockade double junction. Since current through the 
system is accompanied by charging of the grain, interplay between the Coulomb 
forces and the mechanical degrees of freedom can lead to self-oscillations of the 
grain, which in turn, supports charge transport through shuttling of electrons 
between the leads. Much work, both experimental (5, ^ ^ ^ as well as the- 
oretical |i[na[ni[T2[T3[Tl[T3[Ta[niIHlIinillBiniEaE3 has been 
reported in this field. 

In this paper, the noise spectrum of the single-electron shuttle is studied in 
the limit of weak electromechanical coupling. It is found that the onset of me- 
chanical vibrations is accompanied by a peak in the low frequency PSD. Hence, 
by measuring the noise, it can be determined whether or not the grain is oscil- 
lating. This is an important result since direct detection of any high frequency 
mechanical motion is problematic with present day experimental techniques. 

2 Model System 

A schematic picture of the system is shown in fig. QJi. A metallic grain of mass 
M and radius r is placed between two leads, separated by a distance L, via 
elastic, insulating materials. When a bias voltage V = {Vl — Vr) is applied 
between the leads, electron transport occur by sequential, incoherent, tunneling 
between the leads and the grain. The system can also lower its electrostatic 
energy by altering the grain position X . 

The electron transport in the system is described using the notion of the 
equivalent circuit of fig. characterized by resistances and capacitances Rl,r 
and Cl.r which both depend on the grain position X: 

RlAX) = i?o''^exp(±X/A), ClAX) = C^-'^/il ± X/Al,r). 

Here A is the characteristic length scale for tunneling and Cq ' are the capaci- 
tances of the left and right junctions when the grain is at X = 0. The coefficients 



2 



Al,r are typical capacitance length scales for the left and right tunnel junction 
respectively. The grain self-capacitance Co is also included in the model. 

If Ec > h/R{X)C{X), where Ec = e^/2C, tunneUng is well described by 
the "orthodox" theory of Coulomb blockade [53] and the tunneling rates through 
the left (right) junction is 



^ AGfj^{X,V,Q) 



Here (3 is the inverse temperature and AGJ^(X, V,(5) is the decrease of free 
energy as the event {Q, Ql.r) ^ (<3 ± e, Ql,r T e) occurs. The charges Ql,r 
and Q denote the charges accumulated on the left and right leads and the excess 
charge on the grain. The free energy decrease AG J r^^^ ^7 Q) can be expressed 
using the equivalent circuit model of fig. 

Considering the grain motion to be classical and one dimensional we have 
Newton's equation MX = F^xtX-^j^jQ) the grain displacement. The ex- 
ternal force FExt.(-''^, V, Q) acting on the grain includes an elastic force Fk{X), a 
dissipative force F^{X), an electric force F^{X, Q, V) and a vdW force Fydwi^) 
yielding the equation of motion 

MX = FkiX) + F^(X) + Fe(X, Q, V) + F^dwiX). (2) 

For the elastic force we use a phenomenological non-linear potential acting 
as a harmonic well with spring force constant k for small displacements but 
that diverges as a Lennard- Jones potential (12th power) at the positions X^ 
and Xr: 

a b 



where 



k (Xo - Xr){Xo + XlY^ k iXo + XL)iXo - Xr^^^ 

= 



13 Xl+Xr 13 Xl+Xr 

The dissipative force is modelled as a viscous damping term F^{X) = —jX, 
and the electrostatic force is given by 

F,iX,V,Q) - -^Vl + --^Vr - -g — , 

where Q and Vl,r. can be calculated from the equivalent circuit of fig. ^p. Fi- 
nally, the van der Waals force between the grain and the leads is derived from 
a Lennard- Jones interaction potential for the individual atoms |2()| in the leads 
and the grain respectively. For a geometry with a metallic sphere close to a flat 
substrate one finds, 



F,, 



dW 




[k-r + X)\ 
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where i?a,r- are the Hamaker constants for the attractive and repulsive parts 
of the force. Equation ^ together with eq. describe the dynamics of the 
system. 



3 Features of the noise spectrum 

In this section we wiU briefly describe the current noise-spectrum of the model 
system presented in the previous section. Defining the current I through the 
system as the charge transferred per unit time from the grain to the right lead, 
the PSD is given by 



where the brackets denotes ensemble averaging. 27 By doing direct numerical 
integration, i.e. solving the stochastic differential equation eq. (|2J), the current 
as a function of time is obtained and the noise spectrum calculated. A rep- 
resentative result of such a calculation is shown in fig. 13 This spectrum was 
calculated for a system in the shuttle regime. The PSD has been normalized to 
obtain the Fano factor F — Sii/{2e (/)). 

We have divided the spectrum in fig. |21 into four regions marked I-IV. At 
high frequencies, region IV, the Fano factor is close to a value of 1/2 which is 
the value for a static double junction |27]. In region III two strong peaks are 
located at the vibration frequency and the first harmonic. This is a result of the 
periodic charging and decharging of the oscillating grain. The large magnitude 
of the first harmonic stems from the fact that charge exchange between the grain 
and the left lead gives rise to a displacement current in the right lead where we 
measure /. Directly below the peaks, region II, the noise is suppressed below 
the shot noise level of a static double junction, possibly due to the additional 
in time correlations between successive tunnel events induced by the oscillating 
grain. 

The most interesting part of the spectrum however, is the low frequency 
part in region I, where the Fano factor increases. This enhancement of the 
noise is associated with fiuctuations in mechanical energy. In the next section 
we present an analytical treatment of a simplified model system which does not 
take into account non-linear elastic forces, vdW forces or position dependent 
capacitances. Nevertheless, this model suffices to explain the low frequency 
behavior of the noise spectrum in the shuttle regime as well as to predict a 
(quasi) singular behavior of the Fano factor at the instability threshold. 

4 Low frequency noise for weak electromechan- 
ical coupling 

The system is completely described by the conditional probability densities 
p{X, X ,Q,t\Xo, Xo^Qo^to) to find the system with charge Q at time t in the 
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interval [X, X + dX), [X, X + dX) given it was located around (Xq, Xq, Qo) at 
time to- The time evolution of the conditional probability density, as well as 
the unconditional probability density p{X, X, Q, t), is given by the phase space 
equation 



d_ 

di 



p{X,X,Q,t) 



d 



Xp{X,X,Q,t) 



dX 

d ( -fX-kX + F{Q) 
dX \ M 



p{X,X,Q,t) 



rQ,^Qix)p{x,x,Q',t) 

Q'y^Q 

rQ^Q,{x)p{x,x,Q,t). 

Q'y^Q 



(3) 



Here a simplified version of the external force ffext. = — kX + F{Q) where 
F{Q) = Fg{X = 0,Q,V) has been used. Introducing action-angle variables 
{E,(l>) defined through iX,X) = ^2E/Muj'^ (sin(0), u cos{(j))) , where = 
/c/M, eq. ijSll takes the form 



d_ 

dt 



p{E,cj),Q,t) 



_d_ 
^ d(j) 

-{ 

dE 



a; + — sm (p cos ( 

1\T 



-2^Ecos^i 
M 



F(Q)sin0 



2E 

F(g)cos(0) 



p{E,<p,Q,t) 
p{E,cj),Q,t) 



- J2 ^Q^Q'{E,4>)p{E,4>,Q,t). 



(4) 



To facilitate an analytical treatment it is convenient to rewrite eq. Q. We 
thus introduce the dimensionless quantities r = tuj, n = Q/e, fn = F{ne)/F{e), 
£ = El{\Muj^\^), f„^„/ = uj-^Tn^n', e = F{e)/kX, and 77 = 7/(Mwe). The 
probability density can be written in vectorial form as P(£,0, r), where the 
i — i/i component is given by p(5, 0, n(i), t) with ri(i) = (— l)'~^Int [|]. We also 
introduce the matrices 



and, define three operators Ci ~ ~-§i, + f 0) — f , 



sin (f) — Irj sin (f) cos ( 



dtp 

, and £3 = 2^ 



Equation can then be compactly written as 



It]£ cos^ — /Vf cos ( 



|:P(f,<^,r) = 



£1 + 6 £2 + ^3 



p(£:,0,t) 
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Since the low frequency noise is governed by a time scale much longer than the 
period of vibration, adiabatic elimination of fast variables can be used to treat 
slow fluctuations in amplitude. This amounts to doing perturbation theory in 
e. 

Following ref. |2H| we let Pa and Qa be normahzed eigensolutions to £iP\ = 
X{£)P\ and C\Q\ = X*{£)Qx, corresponding to eigenvalue A. Po{£,(j>) is the 
stationary solution to the unperturbed problem, i.e., £iPo = 0, while Qo solves 
the adjoint equation £JQo = 0. The projector Vq onto this state acting on an 
arbitrary vector x(f , 0, t) is given by 

(pox) (£,(/.,r) =Po(£,0)^'"d0'QS(f,0')-x(5,</)',T). 

By direct insertion, using that C\ = + {£,(/)), one finds that (Qo)n = 
[1, 1, • • • , 1]^. Next, splitting up P{£, (f>, t) in mutually orthogonal parts v = 
VqP and w = (1 — Vo)P, an equation, correct to second order in e, can be 
found for v; 

= eVoL^ - (^Vdzt{\L2 + 4 - %U\-w. (5) 

OT 

Here ££4 — I-^ [l{£) — which contains the average dissipation per cycle 

"f{£) and the average pumped energy per cycle yV{£) defined through 

-f{£) = 2eT]£ d<l)Cl*-Po{£,(f>)cos^(t), W{£) = 2eV£ d(j>Q*-fPQ{£,<j>) cosq 
Jo Jo 

Letting p{£, r) be the probability density for finding the system with mechanical 
energy £, i.e. v = Pop, and using eq. ^ a Fokker-Planck equation for p is 
obtained; 

= -^[{Wi£)-j{£)+Oie'))p{£,r)] 

V .fxi£)gxi£) ,^ . 
x^o ^ ' 

where h{^) = \T ^-^Qo ' OP,,{£, ^) and gx{£) = " d^Q^ • OPo(f , 0), with 

6 EE 2t^£ cos2 (j) - cos(/) . Since {yV{£) — 7(f)) is of order e the small 

energy shift of order can be ignored to a first approximation and an equivalent 
Ito stochastic differential equation for the vibrational energy £ is arrived at 

d£ = [W{£) ~ 7(£)] dT + a(£)dVK(r), 




where W(t) is the Wiener process and a(£) = ^\ —2^-^^^ f\{£)g\{£) 
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Performing a small noise expansion [SHI around So, where Sq solves {yV{£o) ~ 
j{£o)) = 0, an analytical expression for the PSD can be obtained 



[ai£o)I'i£o)]' 



Here the current I is given by I{£) = ^ cJ^Qq • J{£ , 4>)Po{£ , (/))■ where 
= '5.jf|^(,)_,„(,)_i(£,0)-f|^(,)_„(,)+i(£: ,0) is the current operator. 

Equation ((BJ is our central result. Two important conclusions can be drawn 
from this equation. First, below a break frequency given by luc — i"f'i£o) — 
W'{£o)) the noise will rise due to fluctuations in oscillation amplitude. This 
agrees with the PSD in fig. |21 Second, approaching the transition point (from 
above) {'y' {£o) — W' {£o)) (see ref. jTU]) which implies that the low frequency 
noise may display a divergent behavior at the point of transition from shuttle 
to stationary regime. 

To investigate this eqs. (|1I2|) were simulated for a sequence of voltages for 
one typical system and the resulting I — V curve as well as the Fano factor 
are shown in fig. |31 The parameter values used are shown in table ^ They 
correspond to a nanometer sized Au grain commonly used in experiments with 
self-assembled Coulomb blockade double junctions. 

Although, as explained in the non-parabolic confining potential smoothens 
any step-structure in the current-voltage characteristics, the transition between 
static- and shuttle-operation is clearly visible in the noise spectrum. Below the 
threshold voltage the Fano factor is of the order 1/2 whereas in the shuttle 
regime, it is higher. In accordance with eq. EJ approaching the threshold from 
above (higher to lower voltages) the PSD shows divergent behavior. 



5 Conclusions 

Using adiabatic elimination of fast variables and numerical integration of eqs. 

the PSD of the classical single electron shuttle has been studied in the case of 
weak electromechanical coupling. We have focussed on the low frequency part 
of the spectrum, which is the part most susceptible to direct measurements, and 
found that the shuttle regime can be distinguished from the stationary regime 
by the noise level at low frequencies. 

In particular, the shuttle regime is associated with a rise in the low frequency 
PSD, as compared to a stationary Coulomb blockade double junction, resulting 
in an increased Fano factor. This increase is due to slow variations in the current 
arising from variations in oscillation amplitude. Approaching the point of the 
transition from above (lowering the bias voltage when in the shuttle regime) the 
noise level shows a (quasi) divergent behavior upon closing in on the transition. 
Hence, even though a measurement of the average current alone may not reveal 
whether the system is in the shuttle regime or not, the accompanying noise 
signature can provide this information. 
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(a) (b) 

Figure 1: Single electron shuttle, (a) A metallic grain of mass M and radius r 
placed between two leads separated by a distance L. The displacement of the 
grain from the center of the system is labelled X. The grain is connected to 
the leads via insulating elastic materials. The leads are biased to the potentials 
Vl and Vr. (b) Equivalent circuit of the system. The tunneling resistances and 
capacitances of the left and right junctions are Rl, Rr, Cl, and Cr. 
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Figure 2: Power spectrum Sniuj) for large amplitude shuttling. For frequencies 
above the vibrational frequency the Fano factor is close to 1/2 as for a static 
Coulomb blockade junction. The peaks are located at the frequency of vibration 
and at the first harmonic. For frequencies below the vibrational frequency, the 
in time correlation due to the periodic grain motion leads to a slight suppression 
of the noise level. At still lower frequencies the the noise is increased due to 
slow fluctuations in oscillation amplitude. 
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V [V] ■ 

Figure 3: Current voltage characteristics plotted together with Sii{lu 0). 
The current is the solid line with the scale on the left ordinate while the Fano 
factor is shown for a discrete set of points with the scale on the right ordinate 
(Lines connecting the points are a guide to the eye). Below the critical voltage 
where there is no sustained grain motion the Fano factor is that of a Coulomb 
blockade double junction. Above the critical voltage the grain is oscillating 
and the Fano factor is increased and shows a divergent behavior at the critical 
voltage in accordance with eq. 10). 



13 



Table 1: Numerical values of parameters used to obtain data of fig. |3| 



Quantity 


Value 


Units 


Quantity 


Value 


Units 


L 


5 


nm 


k 


1 


N/m 


r 


1 


nm 




1 


nm 


M 


10-23 


kg 


Ha 


4.10-19 


Nm 


7 


10-13 


kg/s 


Hr 


10-72 





Quantity Value Units 



A 

pL tdR 

Al,Ar 



0.1 

10 
1 

2.5 



nm 

Mn 

aF 

nm 
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